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Abstract 

Background: Multiple acyl-coenzyme A dehydrogenase deficiency (MADD) is an autosomal recessive disease 
caused by the defects in the mitochondrial electron transfer system and the metabolism of fatty acids. Recently, 
mutations in electron transfer flavoprotein dehydrogenase {ETFDH) gene, encoding electron transfer flavoprotein: 
ubiquinone oxidoreductase (ETF:QO) have been reported to be the major causes of riboflavin-responsive MADD. To 
date, no studies have been performed to explore the functional impact of these mutations or their mechanism of 
disrupting enzyme activity. 

Results: High resolution melting (HRM) analysis and sequencing of the entire ETFDH gene revealed a novel 
mutation (p.Phe128Ser) and the hotspot mutation (p.Ala84Thr) from a patient with MADD. According to the 
predicted 3D structure of ETF:QO, the two mutations are located within the flavin adenine dinucleotide (FAD) 
binding domain; however, the two residues do not have direct interactions with the FAD ligand. Using molecular 
dynamics (MD) simulations and normal mode analysis (NMA), we found that the p.Ala84Thr and p.Phe128Ser 
mutations are most likely to alter the protein structure near the FAD binding site as well as disrupt the stability of 
the FAD binding required for the activation of ETF:QO. Intriguingly, NMA revealed that several reported disease- 
causing mutations in the ETF:QO protein show highly correlated motions with the FAD-binding site. 

Conclusions: Based on the present findings, we conclude that the changes made to the amino acids in ETF:QO 
are likely to influence the FAD-binding stability. 



Background 

Multiple acyl-CoA dehydrogenase deficiency (MADD), 
also known as glutaric aciduria type II, is an autosomal 
recessively inherited disorder of fatty acid metabolism [1]. 
In the majority of cases, the disorder is caused by a defect 
in either the alpha- or beta-subunit of the electron transfer 
flavoprotein {ETFA; OMIM 608053, ETFB; OMIM 
130410) or the electron transfer flavoprotein dehydrogen- 
ase (ETFDH; OMIM 231675) genes. However, in some 
patients, the disorder may result from a few, currently 
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unidentified, disruptions to riboflavin metabolism [2-4]. 
Because of a deficiency of ETF or ETFDH, the acyl-CoA 
dehydrogenases are unable to transfer the electrons gener- 
ated by the dehydrogenation reactions. This results in the 
accumulation of various intramitochondrial acyl-CoA 
esters. Recently, ETFDH mutations were reported to be 
major causes of riboflavin-responsive MADD [5]. 

Characteristics such as organic aciduria may be detected 
only during periods of illness or catabolic stress. In the 
majority of cases, depending on the metabolic status of the 
patient, plasma acylcarnitine analysis may not provide a 
definitive diagnosis. Amongst the most daunting chal- 
lenges facing MADD is metabolic heterogeneity, in parti- 
cular, for the diagnosis of MADD. Thus, numerous 
molecular diagnostic methods have been established for 
analyzing the ETF and/or ETFDH genes for the purpose of 
achieving the final diagnosis of MADD. Previously, we 
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reconfirmed the high prevalence of the C.250G > A (p. 
Ala84Thr) mutation in Taiwanese patients with riboflavin- 
responsive MADD, as well as in the normal population 
using a high-resolution melting (HRM) analysis [6]. It 
should be noted that all the nine patients examined, 
including three pairs of siblings, harbor this mutation and 
four were homozygous. We estimated the allele frequency 
to be 0.004 for the p.Ala84Thr mutant (c.250G > A). 
Thus, the estimated carrier frequency of this mutation in 
the Taiwanese population is 1:125 (0.8%). The homozy- 
gous mutant frequency was 1:62,500. Nevertheless, the 
carrier frequency of C.250G > A is estimated to be 1.35% 
amongst the population in southern China. Thus, 
-1:22,000 Han Chinese are expected to suffer from ribo- 
flavin-responsive MADD [7]. Intriguingly, p.Ala84Thr has 
been identified in Asian countries including Taiwan, Hong 
Kong, Singapore, Thailand and southern China [6-11]. 
A total of 69 patients were examined from the above refer- 
enced literature and diagnosed with late-onset MADD. 
Amongst these patients, the mutation C.250G > A in the 
ETFDH gene was detected in 49 patients with homozy- 
gous mutations and 20 patients with heterozygous muta- 
tions. To date, this mutation has never been reported in 
Western countries [5]. Thus, it has been suggested that 
the most common mutation p.Ala84Thr is ethnic-specific. 
In addition, the mutation of p.Ala84Thr most likely 
accounts for the mild phenotype and response to ribofla- 
vin. Here, a patient presented at 11/2 years represented a 
mild to late-onset MADD and the patient slightly 
improved following riboflavin treatment (Additional file 1, 
Table SI). The patient had compound heterozygous muta- 
tions (p.Ala84Thr/p.Phel28Ser) of the ETFDH gene. Cur- 
rently, there is no related study to explore the underlying 
mechanism for the most common mutation (p.Ala84Thr). 

Flavin adenine dinucleotide (FAD) is a cofactor for the 
electron transfer flavoprotein: ubiquinone oxidoreductase 
(ETF:QO; encoded by the ETFDH gene) which constitu- 
tes the electron-transport pathway for a number of 
mitochondrial flavoprotein dehydrogenases that are 
mainly involved in fatty-acid and amino-acid metabo- 
lism. Riboflavin (the precursor of FAD) treatment has 
been shown to strikingly ameliorate the symptoms and 
metabolic profiles in MADD patients [1,4,6,12]. This 
may enhance the conformational stabilization of the 
mutant ETF:QO protein [13]. Two of the mutations (p. 
Ala84Thr and p.Phrl28Ser) are located in the FAD- 
binding domain; however, the two amino acids do not 
have direct interactions with FAD according to the pre- 
dicted 3D structure of ETF:QO. The closest distances 
between the mutations (p.Ala84Thr and p.Phrl28Ser) 
and the FAD binding residues are 8.6 and 13.6 A, 
respectively. Therefore, to explore the effects of the 
mutations on ETF:QO dynamics, molecular dynamics 
(MD) simulations of the wild type (WT) and mutant 



type (MT) ETF:QO in the same model environment 
were compared. Besides the MD simulations, an alterna- 
tive method, normal mode analysis (NMA), for studying 
protein motions was used to analyze the dynamic corre- 
lations between the mutation sites and the FAD-binding 
motif. 

In this study, molecular modeling, MD simulations 
and NMA were performed to investigate the structural 
implications of the underlying mutations on the protein 
conformation. The obtained results may provide an 
insight into the pathomechanism of MADD. 

Methods 

Patients 

The patient was a girl aged 11/2 years who was clinically 
diagnosed with MADD in 2009. Her dried-blood spot 
sample was sent for acylcarnitine analysis. The results 
demonstrated elevation of acylcarnitines during the illness, 
consistent with glutaric aciduria type 2. Gas chromatogra- 
phy-mass spectrometry revealed significant elevations of 
urinary glutaric acid, 3-methylglutarate, 3-methylcluco- 
nate, adipate, 3-methyladipate, octenedioate, suberate, and 
hexanoylglycine. The present study was approved by the 
Institutional Review Board (IRB) of Kaohsiung Medical 
University Hospital. 

Genetic analysis 

The blood samples were collected from the patient and 
the family members as well as from normal controls 
with informed consents. The genomic DNA samples 
were extracted from peripheral whole blood using the 
NucleoSpin® Blood Kit (Macherey-Nagel), according to 
the manufacturer's instructions. 

Identification of mutations 

HRM analysis was applied for screening the mutation of 
ETFDH, as previously described [6] . PCR was carried out 
in duplicate in a 10 (iL final volume using the LightCycler 
® 480 High-resolution Melting Master (Reference 
04909631001, Roche Diagnostics) 1 x buffer which con- 
tained Taq polymerase, nucleotides, the dye ResoLight 
and 30 ng DNA. The primers and MgCl 2 were used at a 
concentration of 0.25 (iM and 2.5 mM, respectively for 
identifying the ETFDH gene mutations. 

The HRM assays were conducted using the LightCy- 
cler® 480 instrument (Roche Diagnostics) provided with 
the software LightCycler® 480 Gene-Scanning Software 
Version 1.5 (Roche Diagnostics). 

The PCR program required a SYBR Green I filter (533 
nm). The program consisted of an initial denaturation- 
activation step at 95°C for 10 min, followed by a 45-cycle 
program. The 45-cycle program consisted of the following 
steps: denaturation at 95°C for 15 s, annealing at 58°C or 
60°C for 15 s and elongation at 72°C for 15 s with the 
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reading of the fluorescence. The acquisition mode was sin- 
gle. The melting program consisted of three steps: dena- 
turation at 95°C for 1 min, renaturation at 40°C for 1 min, 
and the subsequent melting that included continuous 
fluorescent readings from 60 to 90°C at a rate of 25 acqui- 
sitions per degree centigrade. 

To confirm the results of the HRM analysis, a sequen- 
cing analysis was also performed for all the samples. After 
the HRM analysis, the samples were purified using the 
FavorPrep™ PCR clean-up mini kit. The PCR products 
generated after HRM were sequenced directly. The 
sequence reaction was performed in a final volume of 
10 (iL, comprising 5 [iL of the purified PCR product, 
2.5 \iM of each PCR primer and 1 (iL of the ABI PRISM 
terminator cycle sequencing kit vl.l (Applied Biosystems). 
The sequencing program started from 96°C for 1 min and 
then followed by 25-cycle PCR program. The program 
consisted of the following steps: denaturation at 96°C for 
10 s; annealing at 50°C for 5 s; and elongation at 60°C for 
4 min. The sequence detection was performed using the 
ABI Prism 3130 Genetic Analyzer (Applied Biosystems). 
To ensure that the mutation found in the patient was the 
one specifically present in patients with MADD, 60 healthy 
control chromosomes were also screened. 

Sorting intolerant from tolerant (SIFT) analysis 

SIFT tool [14-16] generates multiple alignments of the 
sequences over different species to look at the conserved 
sequence of a gene; it asses the conserved amino acid posi- 
tions and analyzes the effect of missense changes on the 
conserved structure of proteins over the course of evolu- 
tion. The SIFT tool assigns a score to the mutations, and 
the score of < 0.05 is considered potentially damaging. 

Structure prediction of human wild type (WT) and mutant 
type (MT) ETF:QO 

Three-dimensional (3D) structures of the human ETF:QO 
were predicted using the (PS) server [17,18]. (PS) 2 auto- 
matically uses an effective consensus strategy. The above 
mentioned strategy combines structural- [18] and profile- 
based [19-21] comparison methods for both template 
selection and target-template alignment. (PS) 2 server 
selected the porcine ETF:QO (PDB entry: 2GMH) [22] as 
the model template. The ETF:QO models were built based 
on a sequence and structural alignment from the above 
mentioned template. A total of six models were yielded for 
each target by using a set of parameters. These models 
were based on the alternative target-template alignments. 
The qualities of the modeled structures were checked 
using Ramachandran plots [23]. The model with the high- 
est quality was selected as the final predicted structure for 
the target. Finally, the protein-ligand complex was con- 
structed by superimposing the predicted human ETF:QO 



model to the 3D structure of the porcine ETF:QO-ligand 
complex. 

Molecular dynamics simulations 

Molecular dynamics (MD) simulations were performed 
using the software GROMACS ver. 4.5.1 [24]. OPLS-AA 
all-atom force field was used for the energy calculations. 
The models were solvated with simple point charge 
(SPC) water molecules, and the models were simulated 
in a triclinic box [25] with periodic boundary conditions. 
The simulations were performed in the canonical NVT 
(number of particles, volume and temperature) ensem- 
ble. The energy of the models was first minimized using 
the steepest descent algorithm. This was followed by the 
position-restrained MD simulations for 200 ps. The lin- 
ear constraint solver (LINCS) algorithm was used to 
constrain the bonds [26]. The MD simulations were car- 
ried out at constant pressure and temperature for 2.0 ns 
using an integration time step of 2 fs. The non-bonded 
interactions were cutoff at 10 A. At every 100 time 
steps, the coordinates from the MD simulations were 
saved. 

Analysis of MD trajectories 

The objective was to understand the structural and func- 
tional implications of the amino acid substitution. Thus, 
the trajectories of WT and MT were analyzed for the fol- 
lowing structural properties as a function of time: a) the 
root mean square deviation (RMSD) of the Ca atoms with 
respect to the starting conformation; b) RMSD of the 
FAD-binding motif with respect to the starting conforma- 
tion; and c) the B-factors [27] of the Ca atoms which were 
calculated from the last 1 ns of the MD trajectories. The 
B-factors in protein structures reflect the fluctuation of 
atoms about their average positions. A large B-factor indi- 
cates high flexibility of individual atoms. 

Normal mode analysis 

Normal mode analysis (NMA) is a powerful tool to esti- 
mate the dynamics based on structure [28-31]. Low-fre- 
quency normal modes describing the large-scale real- 
world motions of a protein have been demonstrated to 
be related to biological function [28,32]. In the present 
study, the NMA was performed using the software GRO- 
MACS ver. 4.5.1 by diagonalization of the mass-weighted 
Hessian matrix. The force field-based NMA requires 
careful energy minimization of the initial structure. We 
first minimized the human ETF:QO structure with the 
steepest descent integrator for a couple of steps. Next, 
the L-BFGS integrator was used to minimize the maxi- 
mum force close to zero and subsequently generate a 
Hessian matrix. The NMA was performed using the algo- 
rithm in GROMACS. Finally, the pair-wise residue 



Er et al. BMC Structural Biology 201 1, 1 1:43 
http://www.biomedcentral.eom/1 472-6807/1 1 /43 



Page 4 of 1 1 



correlation was calculated. The above mentioned correla- 
tion was calculated as the average correlation summed 
over 50 lowest nonzero frequency normal modes [33]. 

Results 

Patient and molecular diagnostic 

HRM analysis of ETFDH gene revealed a novel T > C 
point mutation at position 383. The point mutations 
were confirmed by direct DNA sequencing (Figures 1A- 
1D). The hotspot mutation C.250G > A (p.Ala84Thr) 
(Figures 1A and IB) and a novel mutation C.383T > C (p. 
Phel28Ser) (Figures 1C and ID) in the ETFDH gene 
were identified in this patient. These mutations are highly 
conserved amongst the different species (Additional 
file 2, Table S2). The SIFT tool analysis revealed a score 
of < 0.05 and predicted that the replaced amino-acid was 
potentially damaging. This indicates that these mutations 



are most likely pathogenic. Figures IE and IF show the 
melting profiles of the ETFDH mutations in the patient's 
family. The patient had compound heterozygous muta- 
tions (p.Ala84Thr/p.Phel28Ser). The sibling was a carrier 
of C.383T > C (p.Phel28Ser). The father and mother 
were carriers of C.383T > C (p.Phel28Ser) in exon 3 and 
C.250G > A (p.Ala84Thr) in exon 3, respectively (Figures 
IE and IF). The above mentioned mutation was 
not detected in the 60 unrelated control chromosomes 
(wild type) (Figure 1G). The pedigree is presented in 
Figure 1H. 

Structure prediction of human ETF:QO 

Human ETF:QO models were constructed using the 
published porcine ETF:QO crystal structure (PDB Id: 
2GMH) as the template (Figure 2A). A close-up view of 
the hydrophobic residues which are located at the FAD- 
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Figure 1 HRM analysis of ETFDH gene. Electropherograms of various mutations: (A) C.250G > A (wild type); (B) C.250G > A (heterozygous); (C) 
C.383T > C (wild type); (D) C.383T > C (heterozygous); (E) Normalized and temperature-shifted difference plots, the melting profile of C.250G > A; 
(F) Normalized and temperature-shifted difference plots, the melting profile of C.383T > C; (G) The mutation of C.383T > C was not detected in 
60 unrelated control chromosomes (wild type); (H) Affected pedigree with familial segregation of both mutations. Squares, male subjects; circles, 
female subjects. Affected and unaffected subjects are represented by solid and open symbols, respectively. The subject with a sample available 
for the present study is indicated with an asterisk. 
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binding domain, including V71, 173, A84, V85, L87, 
V100, L127, and F128 is shown in Figure 2B. The FAD- 
binding residues are located at the loop that connects 
strand pi and helix al. Therefore, it is reasonable to 




hypothesize that these hydrophobic residues would most 
likely participate in hydrophobic interactions to stabilize 
FAD binding. The ETF:QO mutations of p.Ala84Thr 
and p.Phel28Ser were constructed in the same way as 



A84T F128S 

Figure 2 Structure prediction of ETF:QO. (A) The predicted wild-type model of human ETF:Q0. The structure comprises three domains: the 
FAD-binding domain (gray), the ubiquinone-binding domain (blue) and the 4Fe4S cluster domain (red). The three redox centers: FAD, 
ubiquinone and the 4Fe4S cluster are shown in ball-and-stick representations. A cluster of hydrophobic residues which are located at the FAD- 
binding domain and at a distance below 4 A around the mutation residues (p.Ala84Thr and p.Phe128Ser) are depicted in the sphere model. (B) 
A close-up view of the hydrophobic residues (V71, I73, A84, V85, L87, V100, L127 and F128). The a-helices and p-strands are also labeled. (C) 
Model of the p.Ala84Thr mutant in which the hydrophobic residue A84 is replaced by the hydrophilic residue threonine (yellow). (D) Model of 
the Phe128Ser mutant in which the hydrophobic residue F128 is replaced by the hydrophilic residue serine (yellow). The 3D molecular graphs 
are displayed using PyMOL [51]. 



Er et al. BMC Structural Biology 201 1, 1 1:43 
http://www.biomedcentral.eom/1 472-6807/1 1 /43 



Page 6 of 1 1 



the WT structure. Both replacements are located at the 
center and peripheral of the hydrophobic patch, respec- 
tively (Figures 2C and 2D). 

Molecular dynamics simulations 

For the MD simulations, the trajectories of the WT and 
MT ETF:QO in the explicit solvent were calculated. The 
backbone RMSD values for WT and MT ETF:QO dur- 
ing the production phase relative to the starting struc- 
tures were plotted (Additional file 3, Figure SI). These 
values were plotted to obtain an estimate of the MD tra- 
jectory quality and convergence. For the simulations of 
WT and MT ETF:QO, the obtained data demonstrate 
that after a rapid increase during the first 0.2 ns, the tra- 
jectories are stable with average values of 1.43, 1.38 and 



1.51 A for WT and the p.Ala84Thr and p.Phel28Ser 
ETF:QO mutants, respectively. Statistical analysis of the 
RMSD data reveals that the trajectories are more stable 
after the first 1.0 ns. Therefore, only the second half of 
the trajectory was analyzed. 

Based on Zhang's study [22], a segment corresponding 
to residues 71-90 of human ETF:QO was found to form 
the FAD-binding motif. The binding site was located at 
the pi/al loop (Figure 2B). The RMSD of the FAD- 
binding motif with respect to the starting conformation 
was compared for the WT and MT structures during 
the course of the simulations. The RMSD was found to 
increase as a function of time for the MT p.Ala84Thr 
(Figure 3A) and MT p.Phel28Ser (Figure 3B) when 
compared with the WT. 
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Figure 3 RMSD plots of the FAD-binding motif. Comparison of the RMSD plots of the FAD-binding motif (residues 71-90) of WT and MT (A) 

p.Ala84Thr and (B) p.Phe128Ser structures with respect to the starting conformation during the course of the simulation. WT and MT plots are 
presented in blue and red, respectively. 
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The analysis of the B-factors for each residue at the (31, 
al, a2, and a3 regions revealed that the atomic fluctua- 
tions of the p.Ala84Thr and p.Phel28Ser mutants were 
significant at the pi/al loop (involved in the FAD bind- 
ing) when compared with the WT values (Figure 4). Helix 
a2, which may be involved in the support and stabilization 
of the al and a3 helices, was also found to have higher B- 
factors in both mutants (Figures 4B and 4C). Moreover, 
the fluctuations of the MT p.Phel28Ser were also 
observed to be significant at helix a3 (Figure 4C). From 
the above mentioned analysis, it was observed that the 
mutations of p.Ala84Thr and p.Phel28Ser may induce 
additional fluctuations in the pi/al loop, a2, and a3 
regions. Such fluctuations may promote instability in this 
region of the protein and hamper FAD binding. 

Normal mode analysis 

The correlation map for ETF:QO is shown in Figure 5. 
General structural elements of ETF:QO can be identified 
by the pattern of the cross-correlations. The enlarge- 
ments along the diagonal are typical for a -helices, where 
motions of amino acids adjacent in the sequence are cor- 
related, such as the elements of the helices al, a2 and a3 
(Figure 5). Cross-correlations that extend as plumes from 
the diagonal correspond to anti-parallel elements of the 
secondary structure such as the strands P5, P6, P7 and P8 
(Figure 5). The cross-correlations are also found between 
the multiple structural elements, whereby some result 
from specific contacts, while others are longer range in 
nature. For example, long-range cross-correlations are 
observed between the FAD-binding motif (colored in 



blue) and helix a2 (black square 1). Between the FAD- 
binding motif and helix a2, correlations exist because of 
the presence of a hydrophobic patch consisting of V85 
(helix al), L127 (helix a2) and F128 (helix <x2) (Figure 
2B). Helix a3 also displays significant correlated motion 
with the FAD-binding motif and helix a2 (black squares 
2 and 3). These results suggest correlated dynamics 
between the FAD-binding motif, helix a2 and helix a3. 

Discussion 

In this study, we found that the p.Ala84Thr and p. 
Phel28Ser mutations are most likely to alter the protein 
structure near the FAD binding site as well as disrupt 
the stability of the FAD binding required for the activa- 
tion of ETF:QO using MD simulations and NMA. Inter- 
estingly, NMA revealed that several reported disease- 
causing mutations in the ETF:QO protein show highly 
correlated motions with the FAD-binding site. 

In the FAD-binding family, the pyrophosphate moiety 
binds to the most strongly conserved sequence motif 
[22]. A subset of the proteins that adopt the Rossmann 
fold [34-36] also bind to the nucleotide cofactors such 
as FAD and NAD(P), and function as oxidoreductases. 
The Rossmann fold can often be identified by the short 
amino-acid sequence motif, GxGxxG. The pivotal role 
of the Gly residues in the conserved central GxGxxG is 
well understood [37,38]. The Rossmann fold begins with 
a P -strand connected by a short loop to an a-helix [39]. 
As previously described, the most conserved and well- 
studied sequence motif, xhxhGxGxxGxxxhxxh(x) 8 hxhE 
(D), is part of the Rossmann fold [37,40-42]. Here x is 
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WT A84T F128S 

Figure 4 B-factors analysis. Structures of the (A) WT, (B) p.Ala84Thr and (C) p.Phe128Ser ETF:Q0 proteins are drawn in cartoon putty 
representation at the pi, a1, a2, and a3 regions, where the color is ramped by residue from blue as the lowest B-factor value to red as the 
highest B-factor value. In addition, the size of the tube also reflects the value of the B-factor, where the larger the B-factor the thicker the tube. 
The structures in the other regions are colored in white and displayed in cartoon tube representation, where the size of the tube is independent 
of the B-factors. The red arrows indicate the positions of the mutations. 
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Figure 5 Cross-correlation map calculated from the NMA. A negative value refers to an anti-correlation between residue fluctuations (blue), 
whereas a positive value indicates concerted motion in the same direction (red). 



any residue and h is a hydrophobic residue. The hydro- 
phobic residues provide hydrophobic interactions 
between the a-helix and the P -sheet Intriguingly, the 
mutation of p.Ala84Thr is located at a conserved hydro- 
phobic residue of this motif of the Rossmann fold (Fig- 
ure 2C). In addition, the mutation of p.Phel28Ser is also 
located at this hydrophobic patch (Figure 2D). The 
above mentioned observations, along with the obtained 
results of the trajectory analyses (Figures 3 and 4), 
strongly suggest that these mutations could influence 



the structural stability of the FAD -binding motif. Thus, 
it could lead to unstable FAD binding. Therefore, we 
postulate that p.Ala84Thr and p.Phel28Ser may alter 
the stability of the FAD-binding site and hamper effi- 
cient binding of the ligand. In turn, this instability of 
the FAD-binding site may impair the activity of the elec- 
tron-transport system. 

The binding of FAD is important for the catalytic activ- 
ity of flavoproteins as well as for the correct folding, 
assembly and protein stability [43-45]. A reduced 
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availability of intra-mitochondrial FAD may result in a less 
stable or inactive conformation of the proteins. As a result, 
an accelerated breakdown could result. However, supple- 
mentation of riboflavin is most likely to increase the intra- 
mitochondrial FAD concentration. As a result, FAD bind- 
ing could be promoted. This could ameliorate the effect of 
the mutations that reduce the affinity of ETF:QO for FAD 
[5,13]. It should be noted that riboflavin treatment has 
been shown to ameliorate the symptoms and metabolic 
profiles in recent studies [5,6,8,11,12]. 

If the charge or size of the amino acid is altered and the 
amino acid is buried within the core of the protein, folding 
may not proceed and the result is a complete loss of func- 
tion. Such a loss of function has been found for a large 
number of variant short-chain acyl-CoA dehydrogenase 
(SCAD) and medium-chain acyl-CoA dehydrogenase 
(MCAD) proteins [46,47]. In addition, it may also cause 
conformational changes to the 3D structure by altering its 
hydrophobicity, charge, overpacking, or cavity at the bur- 
ied site, leading to a temperate damage [48]. In the present 
study, we identified a novel mutation, p.Phel28Ser. We 
have used a NMA to characterize the detailed changes to 
protein fluctuations. The results suggest that the dynamics 
of the FAD-binding motif, helix a2 and helix a3 are corre- 
lated. The substitution of p.Phel28Ser may alter the 
hydrophobic interactions between helix al and helix a2 
(Figure 5). It should be noted that numerous mutations 
near residue p.Phel28 have been identified. Wen et al. 
[48] recently identified three ETFDH mutations, p. 
Leul27Arg, p.Aspl30Val, and p.Trpl31Cys in riboflavin- 
responsive MADD patients in the north of China. Er et al. 
[6] and Liang et al. [8] reported one ETFDH mutation, p. 
Leul27His, in riboflavin-responsive MADD patients in the 
Taiwanese population. Law et al. [11] also reported one 
ETFDH mutation, p.Prol37Ser, in a Chinese family with 
riboflavin-responsive MADD. In addition, Goodman et al. 
[49] reported p.L138R as a disease-causing mutation. 
Interestingly, these mutations are located very close to 
each other in the region of the helix a2. Thus, these find- 
ings demonstrate the highly correlated motions of the 
FAD-binding motif (Figure 5 and Additional file 4, Figure 
S2). Previously recognized ETFDH mutations, p. 
Argl75Leu and p.Argl75His [6-8,50], were found in the 
region of helix a3. Both mutations have significantly cor- 
related motions with the FAD-binding motif, as shown in 
Figure 5. Based on the present findings, we believe that 
the mutations to the amino acids in these regions are 
likely to affect the FAD-binding stability. Therefore, it is 
highly recommended that the ETFDH mutation in exon 3, 
4 and 5 be initially screened. However, based on the B-fac- 
tors analysis, the atomic fluctuations of the p.Leul27Arg, 
and p.Argl75His were not significant at the pi/al loop, 
which was involved in FAD binding (Additional file 5, Fig- 
ure S3). According to published literatures, these MADD 



patients are compound heterozygous for two mutations. 
Therefore, it appears probable that other mutations most 
likely account for the disease state. Characterization of 
additional mutations in the other region of ETF:QO will 
be established in future studies. Currently, the present 
study is the first to demonstrate through NMA that sev- 
eral previously reported ETFDH gene mutations have 
highly correlated motions with the FAD-binding site. 

Conclusions 

The present study is the first demonstration of p. 
Ala84Thr and p.Phel28Ser by molecular modeling, MD 
simulations, and correlated motion analysis. Along with 
the obtained results of the MD simulations and NMA, 
the present findings strongly suggest that p.Ala84Thr 
and p.Phel28Ser are likely to alter the ETF:QO protein 
structure near the FAD-binding site as well as disrupt 
the stability of FAD binding, which is essential for the 
activation of ETF:QO. The change of hydrophobic inter- 
action may destabilize FAD binding, giving rise to a less 
stable or inactive conformation. Overall, the obtained 
findings have implications for the pathogenic mechan- 
ism of riboflavin-responsive MADD. 
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